Prenkel-Kontorova Model of Vacancy-Line Interactions on Ga/Si(112): Formalism 
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We describe in greater detail the exactly solvable microscopic model we have developed for ana- 
lyzing the strain-mediated interaction of vacancy lines in a pseudomorphic adsorbate system (Phys. 
Rev. Lett., to appear). The model is applied to Ga/Si(112) by extracting values for the microscopic 
parameters from total-energy calculations. The results, which are in good agreement with exper- 
imental observations, reveal an unexpectedly complex interplay between compressive and tensile 
strain within the mixed Ga-Si surface layer. 

PACS numbers: 68.55.-a, 68.55.Ln, 71.15.Nc, 61.16.Ch 
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When a material is grown pseudomorphically on a 
lattice-mismatched substrate, the resulting strain field 
can lead to self-organized structures with a length scale 
many times the atomic spacing. One well known exam- 
ple is the Ge/Si(001) dimerized overlayer system. The 
Ge film is compressively strained (by 4% relative to the 
bulk), and the system lowers its energy by creating dimer 
vacancies in the surface layer; at the vacancy sites, the 
exposed atoms in the second layer rebond to eliminate 
their dangling bonds. The missing-dimer vacancies order 
into vacancy lines (VLs) with 2 x TV periodicity, where 
the optimal N depends on the Ge coverage @. Even 
for coverages as low as a few monolayers, the concept 
of elastic strain relaxation within a coherent pseudomor- 
phic Ge film is appropriate. For example, Tersoff showed 
theoretically for a 3-layer Ge film that the equilibrium 
N corresponds to the vacancy density at which the com- 
pressive stress from the Ge overlayer cancels the tensile 
stress from the rebonded missing dimers S . 

For monolayer and lower coverages, the concept of 
strain relaxation becomes problematic, because the strain 
within a partial overlayer becomes difficult to define. In 
this Letter we develop a model for analyzing such situa- 
tions and apply it to another VL system, Ga on Si(112). 
We show that despite striking similarities in the phe- 
nomenology, the underlying energetics of Ga/Si(112) is 
rather distinct from Ge/Si(001). To make our treat- 
ment physically transparent but also quantitatively ac- 
curate, we develop an exactly solvable model of VL in- 
teractions in which the microscopic parameters are ex- 
tracted directly from first-principles total-energy calcu- 
lations. This model contains only nearest-neighbor har- 
monic interactions but reproduces the first-principles re- 
sults quite accurately, and thus allows for a particularly 
simple analysis of the dominant interactions. 

When Ga is deposited on Si(112) and annealed, a 
well-ordered surface is formed consisting of large (112)- 
oriented domains, as shown in Fig. |j} The VLs are ori- 
ented horizontally in Fig. 0. Even at room temperature, 




FIG. 1. Empty-state STM image (300 x 300 A) of 
Ga/Si(112). The wedges mark the vacancy lines, which are 
oriented along the [111] direction (and are interrupted by a 
narrow monolayer-high terrace in the center of the figure) . 



the VLs show minimal thermal meandering — only single 
kinks (up or down by one lattice spacing) are observed 
in scanning tunneling microscopy (STM). The mean VL 
spacing in Fig. |l| is 5.2, with the distribution sharply 
peaked around 5 and 6. The microscopic structure of VLs 
on Ga/Si(112) is shown in Fig. ||. Since Ga is trivalcnt 
it prefers to adsorb at three-fold surface sites. The bulk- 
terminated Si(l 12) substrate, which may be regarded as 
a sequence of double- width (lll)-like terraces and single 
(lll)-like steps, offers just such three-fold sites at the 
step edges, as shown in Fig. ||(b). A single Ga vacancy 
leaves two step-edge Si atoms exposed, which rebond to 
form a dimer. This model was first proposed by Jung et 
al. |H and subsequently confirmed by Baski et al. Q us- 
ing STM and total-energy calculations. It predicts that 
two adjacent Ga vacancies will be very unlikely — because 
the three Si atoms exposed can form only one dimer, re- 
sulting in an extra Si dangling bond; indeed, there are 
no adjacent vacancies visible in Fig. In this way, VLs 
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FIG. 2. (a) Side and (b) top views of Ga/Si(112) with va- 
cancy period N=b. The fully relaxed coordinates are from 
first-principles total-energy minimization, (c) Bonding chain 
of Ga (dark) and Si (light) atoms, (d) One-dimensional 
Frenkel-Kontorova model representing this Ga-Si chain. The 
(horizontal) atomic displacements, Uj, are defined relative to 
the ideal positions (shown by thin grid lines) of the substrate 
atoms; vertical displacements represent the individual sub- 
strate-strain energies, (l/2)fcju|, appearing in Eq. (|l|). 



on Ga/Si(112) are essentially constrained to have a fixed 
width of one vacancy. 

To construct a microscopic model for the interactions 
between VLs we make three simplifying assumptions 
about their ground-state structure. (1) We take the VLs 
to be perfectly straight; this turns a two-dimensional sur- 
face problem into an effective one-dimensional system. 
(2) We assume this one-dimensional system to be periodic 
with a single vacancy per unit cell, that is, the vacancy 
separation is taken to be L = Na, where N is an inte- 
ger and a is the surface lattice constant. (3) Full struc- 
tural relaxation within the local-density approximation, 
described below, shows that the Si substrate atoms are 
essentially unperturbed from their ideal locations (Fig. || 
shows the relaxed geometry for N = 5; substrate atoms 
are white). Thus, we consider the vacancy- vacancy in- 
teraction to be mediated entirely by the bonding chain 
of Si and Ga atoms shown in Fig. ||(c). 

We now map this Ga-Si chain (with vacancy period N) 
onto a one-dimensional chain of harmonic springs con- 
necting N Si atoms and N — 1 Ga atoms in a unit cell of 
length Na, as shown in Fig. ||(d). The bonding of each 
atom, j, to the substrate is also taken to be harmonic, 
with potential minima at the ideal substrate positions. 
This is a variant of the Frenkel-Kontorova (FK) model, 
which has been widely used to study complex static and 
dynamical phenomena arising from purely local interac- 
tions. Within the FK model, the total (potential) energy 
is 
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Here the one-dimensional atomic displacements Uj are 
defined with respect to the ideal positions of the substrate 
atoms, (a/2)j, and the corresponding spring lengths bi 
are defined in Fig. ||(d). This FK model has six parame- 
ters: within each unit cell, one spring represents the Si-Si 
bond at the vacancy (with spring constant Ksi and equi- 
librium length fog?), while the remaining 2N — 2 springs 
represent the Ga-Si bonds (with spring constant K and 
equilibrium length 6 cq ). Two different substrate spring 
constants, fcsi and fcca, represent the bonds from Si and 
Ga atoms in the chain to the rigid substrate. 

To solve for the displacement field that minimizes the 
FK energy, we first consider the force equations for atoms 
away from the vacancies, and then apply the boundary 
conditions due to the vacancies. In equilibrium, the force 
on atom j is given by 

Fj = = K(2v,j — Uj-i — Uj+i) + kjUj, (2) 

which can be written for the Si or Ga sublattice as 

[{2K + k Ga )(2K+k si ) ~ 2k 2 ] Uj = K 2 (u. J+2 +u^ 2 ) . 

(3) 

In the continuum limit (large N) this becomes a simple 
differential equation, 



[f3-l]u(x) = (a/2) 2 u"(x), 



where we have defined 
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Equation (Q) shows that the continuum displacement 
field, u(x), has an exponential solution with the decay 
length (a/2)/y/P — 1. Note that for the case of weak 
binding to the substrate, kj <C K, this decay length 
reduces simply to yj K/k (a/2), where k is the average 
substrate potential. 

For the discrete case (arbitrary N) , Eq. (||) provides a 
recurrence relation whose solutions have the form 



_ Si,Ga 2\j _|_ Si.Ga — 2Aj 

where e ±2A are given by 



e ±2A = yfj>±yfj- 



1. 



(6) 



(7) 



The general form of the discrete solutions is again seen to 
be exponential, although with a somewhat more compli- 
cated form for the (dimensionless) decay length, (2A) _1 . 
It is easy to verify by Taylor expansion that for weak 
substrate binding we again recover the correct limiting 
behavior, (2A)" 1 -> y^Kjk. 
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Eq. (g) provides the following relationship between the 
four coefficients c^ 1 and c„, which are used for the Si and 
Ga site respectively: 



' 1 + k Si /2K 



(8) 



It is convenient to introduce a dimensionless parameter, 
A, defined by equating the rhs of Eq. @ with exp(2A). 
The deviation of A from zero measures the relative differ- 
ence in strength between Si- and Ga-substrate binding, 
and will give rise (see below) to oscillations in the strain 
field around the simple exponential behavior found in the 
continuum limit. 

The remaining unknowns are determined by the 
boundary conditions at the vacancies. By applying these 
we arrive at the exact closed-form solution for the dis- 
placement field and corresponding bond lengths. For ex- 
ample the Si-Si vacancy-bond length, as a function of 
vacancy period, can be written as 



bo(N) = a - 



2uf 



1 + £{coth[2(AT- 1)A] - 1} ' 



(9) 



where uf and £ (both positive) are combinations of the 
various FK parameters. The Ga-Si bond lengths have a 
more complicated form due to the two different substrate 
potentials: 

bi(N) = - + 4c {cosh A sinhA cosh[(2i - 2N + 1)A] 

± sinhA cosh A sinh[(2i - 2N + 1)A]}, (10) 

where the last term (present only when A ^ 0) is pos- 
itive and negative for even- and odd-numbered bonds, 
respectively. The magnitude of the strain field is given 
by the prefactor c = (l/2)(6 - a) e~ A csch [2(1 - AT) A]. 
It is evident from Eqs. (||) and (|l(]) that the strain field 
is characterized by a single length scale, (2A) _1 , which 
describes both the relaxation of each bond with respect 
to the distance, i, from a vacancy, and the relaxation of 
all strains with respect to the vacancy period, N. 

To apply this general solution to Ga/Si(112) we must 
determine numerical values for the six FK parameters 
(two bond lengths, b eq and bf^, and four spring constants, 
K, K$i, fcca ; kgi) appearing in Eq. (|l|). We do this either 
analytically from the Stillinger- Weber potential (for Si- 
Si parameters) or numerically from first-principles total- 
energy calculations (for Ga-Si parameters). These calcu- 
lations were performed in a double-sided slab geometry 
with six layers of Si and a vacuum region equivalent to 
five layers of Si. Total energies and forces were calculated 
within the local-density approximation (LDA) with gra- 
dient corrections using Troullier-Martins pseudopo- 
tentials and a plane- wave basis with a kinetic-energy cut- 
off of 8 Ry, as implemented in the fhi96md code M . To- 
tal energies were completely converged with respect to 
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FIG. 3. Equilibrium bond lengths from the FK model and 
LDA. Left: Si-Si vacancy bond length vs. vacancy period. 
Right: Ga-Si bond length vs. "bond number" [relative dis- 
tance from the vacancy bond, as defined in Fig. 0(d)]. 

Brillouin-zone sampling. Full structural relaxation was 
performed on all atoms except those in the innermost 
double layer until the surface energies were converged to 
0.1 mcV/A 2 . 

Using the equilibrium structure of Ga/Si(112) with no 
vacancies we find that, for small displacements, fcgi = 2.0 
eV/A 2 and k G& — —0.7 eV/A 2 . The negative spring con- 
stant here indicates that in the absence of defects the Ga 
sublattice is at a point of unstable equilibrium — a finding 
confirmed by our LDA results for finite N (below). From 
calculations on isolated infinite Ga-Si chains, we obtain 
K = 9.8 eV/A 2 and b cq = 2.00 A. This equilibrium bond 
length implies a 4% compressive epitaxial strain with re- 
spect to the Si substrate, and contributes to the driving 
force for vacancy formation in Ga/Si(112). Finally, by 
expanding the radial part of the Stillinger- Weber poten- 
tial about the LDA vacancy bond length for N = 2, we 
obtain K Si = 5.5 eV/A 2 and tijg = 2.25 A. 

We now compare the structural predictions of the FK 
model to those of first-principles calculations. Using the 
LDA approach described above, supercells with vacancy 
periodicities from N = 2 to 8 were completely relaxed. 
Fig. H shows the FK and LDA results for the vacancy 
bond length vs. N and for the Ga-Si chain bond lengths 
for N — 8. The agreement is remarkably good, especially 
considering the simplicity of the model. Note the oscilla- 
tory behavior of the Ga-Si FK bond lengths — due to the 
unequal Ga- and Si-substrate binding strengths — which 
is strikingly confirmed by the LDA results. 

We turn next to the description of strain energetics 
and, thereby, the effective VL interaction. Within the 
FK model, we use Eqs. @ and ( ^p| ) to evaluate the strain 
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FIG. 4. Relative surface energies, 7(iV), vs. vacancy den- 
sity, calculated within the FK model and LDA. 

energy, U, for arbitrary N. By exploiting the fact that 
the sums over atomic sites are geometric series, we obtain 
a closed-form solution for U(N) (see Appendix). It is 
important to note that energies for different N are not 
directly comparable, since the number of Ga atoms per 
unit length varies with N. In general, this is resolved by 
considering the surface free energy (per unit area), 

j(N) = (NA)- 1 [E t (N) — (N — l)/x Ga ] /2, (11) 

where E t (N) is the (LDA) total energy for a Ga/Si(112) 
cell with vacancy period N, and /ZQa is the Ga chemical 
potential (the energy per atom of bulk Ga). The chemical 
potential is not derivable within the FK model, but must 
be included "by hand." We do this by defining 7 FK (7V) 
analogously to Eq. ( [TT|) , with a fictitious chemical poten- 
tial adjusted to give identical VL formation energies, Sf, 
within the FK model and LDA (which gives £/ = —43 
meV/A). 

In Fig. [| we show the resulting FK and LDA surface 
energies, relative to the limit of large N. The LDA energy 
minimum occurs at N = 6, in excellent agreement with 
the experimental results for this system, while the FK 
minimum is only slightly higher, at N — 8. The energy 
scale for the effective VL interaction is extremely small: 
relative to infinite separation between vacancies, the LDA 
energy minimum is only of order 2 meV/A 2 , while the FK 
result is about half this value. In general, the agreement 
between the two results is quite remarkable, considering 
the length scale of the vacancy spacing, the smallness of 
the energy scale, and the extreme simplicity of the FK 
model. 

By analyzing the individual contributions to the FK 
strain energy, it is now simple to identify the domi- 
nant physical mechanism determining the equilibrium VL 
spacing. In Fig. |J we plot (as small dots) only those 
terms in Eq. (yj) representing strains within the Ga-Si 
bonding chain, and exclude the contribution from the Si- 



Si rebonded dimer as well as all interactions between the 
chain and substrate atoms. It is clear that this Ga-Si 
chain strain completely dominates the energetics for all 
physically relevant densities. This result has two impor- 
tant implications. (1) The role of the rebonded Si dimers, 
although obviously crucial for eliminating the dangling 
bonds created by the Ga vacancies, plays no further sig- 
nificant energetic role in determining the equilibrium VL 
spacing in Ga/Si(112). This is quite different from the 
role of the rebonded Ge-dimer vacancies in Ge/Si(001), 
for which the tensile stress contribution has been shown 
to cancel the compressive contribution from the overlayer 
only at the proper density ||. (2) While a low density 
of Ga vacancies allows for relief of compressive strain, 
when their density becomes too high, part of the Si-Ga 
bonding chain experiences tensile strain — which acts as 
a repulsive interaction between VLs. This effect is read- 
ily visible in the right panel of Fig. ||. Relative to the 
equilibrium Si-Ga bond length of 2.00 A, bonds far from 
the vacancy (numbered 5-7) are compressively strained, 
bonds near the vacancy (3-4) are essentially unstrained, 
but bonds very close to the vacancy (1-2) are tensilely 
strained. This tensile contribution begins to dominate 
the surface energy at about N < 5. 

In summary, we have developed a Frenkel-Kontorova 
model to analyze the microscopic origins of vacancy-line 
interactions on Ga/Si(112). The model reveals that the 
mechanism of strain relaxation in this submonolayer ad- 
sorbate system is quite complex. In particular, we have 
identified the microscopic origins of attractive and repul- 
sive interactions between vacancy lines: both are medi- 
ated by a combination of compressive and tensile bond 
strains within a single chain of Ga and Si atoms. The sum 
of these strain energies is minimized at a vacancy-line 
density very close to the experimentally observed value. 
In general, we expect the future analysis of similar strain- 
induced self-organized adsorbate systems to benefit from 
this type of simple but accurate analytical model. 

Computational work was supported by a grant of 
HPC time from the DoD Major Shared Resource Cen- 
ter ASCWP. This work was funded by ONR. 

APPENDIX: SOLUTION OF THE 
FRENKEL-KONTOROVA MODEL 

1. The Model 

Let the atoms be located at Xj for j = 1, . . . , 2N — 1. 
Define the displacement by uj = Xj — (a/2)j. Then the 
potential energy is 

2N-2 2JV-1 

U= E ^Ki(h-b?) 2 + £ \k,ul (Al) 

t=0 j=l 
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where Ki — K and b cq = b eq for i ^ and Kq = Kgi and b^ q = &si- The bond lengths are defined as 

2 



X i+ i - Xi 
Xl — X-i 



u l+1 - Ui + § 
ui — U-i + a 



fori = l,2,...,2iV- 
for { = 0. 



In the bulk (away from the vacancy) the force experienced by an atom in equilibrium is 

Fj = = K(2uj — Uj—i — Uj+i) + kjUj 

= K 



k ■ 

2 + ^)Uj- Uj-l - u 



3+1 



At the vacancy, the force on the first Si atom is 



Fi = = K Si { Ul -u-i + a- bsi) ~ K(u 2 - m + - - b cq ) + k Si u v 
Assuming periodicity, the chain is reflection symmetric about j = and j — N, so that 

Uj = U2N—j 

Ul = -U2JV-1 = -U-l. 

Clearly we need only consider one end of the chain. ( |A5| ) then reduces to 

= (2K Si +K + ksi)u 1 -Ku 2 + (2K Si - K)^ + {K b cq - K Si b Si ) . 



(A2) 

(A3) 
(A4) 

(A5) 

(A6) 
(A7) 

(A8) 



2. The Displacements 



The bulk force equation (A4) is linear, so the solution is a linear combination of exponential terms such as 

Uj = c 3 a 3 . 



(A9) 



We plug this into (A4) with the relation Cj — Cj_2, which says that the Ga and Si displacements are scaled differently, 
but that all of the displacements for cither type of atom arc scaled by the same factor, that is, Cj is csi for odd j and 
CQ a for even j . The result is 



2 + -g I Ci = I a + - ) c 



(A10) 



First, we solve for a. Dividing the equation (A10) on two adjacent sites and using Cj — Cj-2 to eliminate a, we get 

(All) 



fcsi 
2K 



2K 



where kj is fcsi for odd j and fcca for even j. This may be solved for a: 



(3= 1 



fcsi 
2K 



1 



2K 



(A12) 
(A13) 



Note that a + = \jo.-. Also note that a could be complex for negative kj. Basically this says that an inverted 
substrate potential can in principle cause the system to become unstable. If we assume kj <^ K, which is true for our 
system, then the system is stable as long as fcsi + &Ga > 0. 

Next, we solve for Cj/cj-!. Multiplying the equation (A10) on two adjacent sites and using Cj — Cj_2 to eliminate 
the Cj's, we get 



CGa 
CSi 



' 1 + k Si /2K 
l + k Ga /2K' 



Thus, the displacements are given by 



p-V 2 c+ ai~ N +p- 1 /2 c _ a N-j 
pV 2 c+ ai~ N + p 1 ' 2 c_ a N ~i 



for j odd (Si), 
for j even (Ga) 



(A14) 



(A15) 



with a = a + in (A13) and p given by (A14). 
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3. The Vacancy Bond Length 



This leaves two undetermined coefficients, c+ and c_, which are set by the boundary conditions at the vacancies. 
As noted above, the periodicity induces a reflection symmetry (A£), so that 



-c_ = c 



(A16) 



and there is only one undetermined coefficient. Substituting the bulk displacements ( A15 ) into the vacancy force 
equation (|A8|) we get 



(K Sl -\K)a + {Kb c *-K S{ b si ) 



' -p- 1 / 2 (2K Si + K + k si )(a 1 ~ N - a N - r ) +Kp 1 / 2 (a 2 - N - a N ~ 2 ) ' 
Then the bond length at the vacancy is given by 

(2K Si - K)a + 2{K - K si b Si ) 



bg = 2ui + a = a + 



-(2K Si + K + k Si ) + Kp ( gi^I^Ii ) ' 



(A17) 



(A18) 



4. The Bond Lengths 

We can rewrite the displacements ( A15| ) in the form 

Uj = l c | [(pV2 + p -l/2) + (p l/2 _ p -l/2 } cog(7ri) j (a J-N _ a N-j^ 

The bulk bond lengths are computed from (A15), ( A16| ) and (A17): 



u i+ i — m + o< 



a c 
2 + 2 



p 1/2 +p- 1/2 ) (a^ 2 -n. 



a i-N+l/2 + a N-i-l/2 



" 2 ) ( 

(t«) (p 1/2 - P - 1/2 ) (a 1 ' 2 + a' 1 / 2 ) (a i - JV+1 



= - + 4c [cosh A sinh A cosh([2i -2N + 1]A) 



cos(7n) sinh A cosh A sinh([2i - 2 A + 1]A)] 



where A = \ log a and A = \ log p. 



5. The Energy 



(A19) 

(A20) 
(A21) 

(A22) 



Now that we have the expressions for the displacements (A19) and the bond lengths (A22), they can be plugged 
into (Al) to compute the energy, U . The following identity is useful to simplify the geometric series that occur in the 
expression for the energy: 



JY 



-N 



x N+ 2 — x N i 



I! — X 2 



Also it is useful to reexpress the displacements and couplings: 



|c { [(pV 2 + p-^ 2 ) + ( P V 2 - p- l ' 2 ){-iy] (a*- N - a""')} 



2 12 



[{p + p- 1 ) + (p- p- l )(-iy] (a 2 >- 2N + a 2N - 2 > 2) 



\ [(fcG a + fc S i) + (fc G a-A;si)(-l) j ] 



(A23) 



(A24) 

(A25) 
(A26) 



Using (A25) and (A26), an expression for the substrate bond energy may be derived by segregating the terms into 
those that alternate in sign, and those that do not. Then the identity (A2S) may be used to simplify the series: 



G 



2JV-1 
3=1 



-(2JV-1) + 



-1 



1 



-C 2 (fcGaP + feiP 



1 



-(2JV-1J + 



a 2N-l +a l-2 Nl 



a + a 1 
sinh (2A[2iV - 1]) 



^c 2 (fc Ga p - fcsiP X ) 



2c 2 if (cosh 2A - cosh 2A) 



sinh (2A) 
1 cosh(2A[27V- 1]) 



+ 2c 2 K sinh(2A) 



-1 



cosh (2A) 

[ ' sinh(2A) 

cosh(2A[2iV- 1])" 



cosh (2A) 



Now consider the bond energy. Let 



ao = \a - 6 oq , 

C\ = 4c cosh A sinh A, 

C2 = 4c sinh A cosh A. 

Then the bond stretch (and its square) away from the vacancy are 

Shi =bi- 6 cq = a + ci cosh ([2i - 2N + 1]A) 
- (-1)^2 sinh ([2i -2N + 1]A) 

Shi = eg + c 2 { | + | cosh(2[2i - 2iV + 1]A)} 
+ c 2 {-i + ±cosh(2[2i - 2JV + 1]A)} 

- c 2 (-l) 4 {ci sinh (2[2i -2N + 1]A) + 2a sinh ([2i - 2N + 1]A)} 
+ 2a oCl cosh ([2i - 2iV + 1]A) . 



The following identities are useful: 



2JY 



^cosh(K [n-N-i]) 

71=1 

^(-l) n smh(i? [n- JV- §]) 



2iY 



s>nh(NK) 
smh(^K) 

sinh(NK) 
cosh(iiT) ' 



They are essentially geometric series like ( A23 ). Also we use the following expressions 



c\ + cj = 8c 2 [cosh(2A) cosh(2A) - 1] , 
c\-c\ = 8c 2 [cosh(2A) - cosh(2A)] , 
Cl c 2 = 4c 2 sinh(2A) sinh(2A), 
cosh(A) - sinh(A) = e~ A . 



From (All), we have 



cosh(2A) 



fcsi 

2a: 



2a: 



(A27) 



(A28) 



(A29) 



(A30) 
(A31) 
(A32) 



(A33) 



(A34) 

(A35) 
(A36) 



(A37) 
(A38) 
(A39) 
(A40) 



(A41) 



The non- vacancy bond energy is then computed from (A34) 
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2N-2 f 

J2 -K8bl = \K\\{2al + cl-cl)(2N-2)- 
»=i ^ 



^ a , ^sinh(4A[A-l]) 



K c i + C 2) 



c\c 2 



2aoCi 



sinh (2A) 
sinh (4A[JV- 1]) 

cosh (2A) 
sinh(2A[iV- 1]) 



2a C2 



sinh (2X[N - 1]) 
cosh (A) 



sinh (A) 

\K 6 cq ] 2 + 4c 2 [cosh(2A) - cosh(2A)])(2A - 2) 

sinh (4A[iV- 1]) 



4c 2 [cosh(2A) cosh(2A) - 1] 



sinh (2A) 



4c 2 tanh (2A) sinh (2A) sinh (4A[A - 1]) 
(|a - 6 cq )ce- A sinh (2\[N - 1]) J . 



The vacancy bond energy is given by 



±K Si 5b 2 = i^si { (a - bsi) 



(2K Si -K)a + 2(Kb^-K si b Si ) 
< 2K Si + K + ksO+Kp( ^\^ ) 



Finally, 



(Ksi- ±K)a+(Kb°*-K Si bsi) 



2p- 1 / 2 (2K Si + K + fc S i)sinh(2A[AT-l]) - 2K p 1 / 2 sinh(2A[iV-2]) 
\ p" 1 ' 2 { [{K si /K) -\]a+ [V* - (K Si /K) hi] } 
p- 1 {2{K Si /K) - l)sinh(2A[A-l]) + sinh(2AA)' 



where we have used 



fc Ga = 2K (p- 1 cosh(2A) - 1) 
k si = 2K (p cosh(2A) - 1) . 



The energy, U, is the sum of ( A29 ), (A43 ) and (|A44| ), using (|A46|) for c. 



6. Simplifications 



Below is a collection of simplified formulas. We can rewrite (A18) in the form of Eq. 

2uf 



b = a 



where 



and 



1 + £ {coth[2(A - 1)A] - 1}' 

(K Si ~ K/2)a + (K b c * - K Si b Si ) 
(2K Si + K + k Si )-Kpa-i 

Kp sinh(2A) 



* -(2K Si + K + k Si ) + K pa- 1 ' 
Further, we can express the coefficient c in terms of bo, rather than vice- versa: 

c = |Oo - a) e- A csch (2[1 - AT] A) . 



* On leave at: Fritz-Haber-Institut, Faradayweg 4-6, D-14195 Berlin-Dahlem, Germany. Electronic address: 

erwin@dave.nrl.navy.mil 
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